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I. INTRODUCTION 

The main purpose of the line of work outlined in this 
conference, is to establish the relationship between grav¬ 
itational radiation and source properties. Thus, for ex¬ 
ample, we known that gravitational radiation is an irre¬ 
versible process, accordingly there must exist an entropy 
production factor in the equation of state (dissipation) 
of the source. 

Since we are dealing with gravitational radiation, we 
need to depart from the spherical symmetry. On the 
other hand, we shall rule out cylindrical symmetry on 
physical grounds. Thus we are left with axial and re¬ 
flection symmetry, which as shown in |3 is the highest 
degree of symmetry of the Bondi metric [6| , which do not 
prevent the emission of gravitational radiation. 

We are using the 1-1-3 formalism 01 , in a given 
coordinate system, and we are going to ressort to a set of 
scalar functions known as Structure Scalars 0 , which 
have been shown to be very useful in the description of 
self-gravitating systems [lll - l^ . 

II. BASIC EQUATIONS, CONVENTIONS AND 
NOTATION 

We shall consider fluid distributions endowed with ax¬ 
ial and reflection symmetry, and we shall assume the line 
element to be of the form: 

= -A^df+B"^ {dr"^+ r'^de'^)+C^d(j)'^ + 2Gdedt, ( 1 ) 

where A,B,C,G are positive functions of t, r and 9, 
and coordinates are numbered as: = t,x^ = r, x'^ = 

9, x^ = 4>. 


The energy momentum tensor describes a dissipative 
fluid distribution and in its canonical form may be writ¬ 
ten as: 

Tap = (M + P)yaVp + P <?q/ 3 + Ila/S + QaVp + QpVa- (2) 
with 

U = TapV^VP, qa = -piVa - TapY^, (3) 


p = \h‘^^Tap, liap = (T^, - , (4) 


hfj.u = gujy + Yi^Vfj,, (5) 


U“ = (2,0,0,0); Va = {-A,0,j,0). (6) 

where pi, P,Ilap,qa,Va denote the energy density , the 
isotropic pressure, the anisotropic tensor, the dissipative 
flux and the four velocity respectively. 

Next,in order to form an orthogonal tetrad, let us in¬ 
troduce the unit, spacelike vectors K,L, S, with compo¬ 
nents 

42 n2„2 I r<2 

A„ = (0,R,0,0); L„ = (0,0,^^- ^ ,0), (7) 


,S„ = (0,0,0,C), (8) 

satisfying the following relations: 

VaV^ = -K°^Ka = -L°^La = -S^^Sa = -1, (9) 

= V°^La = V°‘Sa = 

K°‘La = K°‘Sa = S°‘La = 0 . 
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( 10 ) 
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In terms of the above vectors, the anisotropic tensor may 
be written as 

= i(2n/ + nii){K^Kp - 


where 


with 


+ l{2Iln + nj){L^Lp-^) 

+ ‘^^KLK(^aLi3), 

Ukl = , 


( 11 ) 

( 12 ) 


Hi = {2K°‘KP - - S°^Sf^)Tc,p, (13) 

Hu = {2L°‘Ll^ - - K°‘Kl^)T^p. (14) 

For the heat flux vector we may write 

9^ = qiKf^ + qiiL^, (15) 

or 

qiiG qi Aqii 


q^ = i 


AVA^B^r^ + G2 ’ B ’ ^A^B^r^ + 


,0), (16) 


/ + G^qii \ , , 

qi.= i0,Bqi,A ---^,0 . (17) 


A. Kinematical variables 

The kinematical variables (the four acceleration, the 
expansion, the shear tensor and the vorticity) are defined 
respectively as: 


= V^Va-p = ajKa + aiiLo 


= 0 


e = v.‘i 


An ^ 

^’(42 




I G^t 


+ ) ’ (^®) 




^ 2^252 Q2 

^ G2 fB,t A,t . G.t . G. 


A^B^ \ B A G G 


(19) 


O'a/3 — ^a;/3) “1“ (^0) 


or 


^a/3 — 3^^^^ j3 

+ ^(2(77/ + (Jj){LaL0 — —hajs): 


( 21 ) 


9 ^ 3/5, C, 

2a, + a,, = -(—- — 


( 22 ) 


2(77/ + CT/ = 


AB\-^ 

B G 


A2B2r2 + G2 

^ Gt _£t 

A I A G G 


(23) 


Wr 


= ^ V'^ = ^ 77a/3;.. 1"", (24) 


where fia /3 = V[a-p] + 0[aVs], Wa and rja/Situ denote the 
vorticity tensor, the vorticity vector and the Levi-Civita 
tensor, respectively; 


lla/3 — 


CUqi - ^\Sqj . 


H = 


^( C?,r \ 

G A ) 

2BVA2B2r2 +G2' 


(25) 


(26) 


(27) 


Observe that from (EH) and regularity conditions at 
the centre, it follows that: G = 0 <tt> fl = 0. 


B. The orthogonal splitting of the Riemann Tensor 
and structure scalars 


Using the well kown decomposition of the Riemann 
tensor in terms of the Weyl tensor, the Ricci tensor and 
the Ricci scalar, and linking the two later variables with 
the energy momentum tensor, via the Einstein equations, 
it can be shown that the Riemann tensor may be written 
as: 

^ i^S — ^(F) ^(Q) ,yS AE) ,yS -^(H) lyS’ 

with 

R 


- —(m + 3P)U[“l^./7f; + /7fj,(29) 


{F) v5 2 


Rt) = -167rU[“/ifJg5] - HSFV^.h^^q^^ - (30) 


-16^U[“U[,nfJ + 167r/i[“nfJ 
<) = 4U[“U[.ufl + 4/7[eEfl, (31) 

RIh) .8 = (32) 
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In the above, i?a/3, Hap denote the electric and magnetic 
parts of the Weyl tensor, respectively, defined as usual 
by: 

( 33 ) 

where eapp = ’n^aPpV’'. 

In our case these tensors may be written in terms of 
five scalar functions as: 

EaP = 2 ( 2 ^/ + £ll){KaKp — -hap) 

+ 2 ( 2 ^// + £l){LaLp — —hap) 

+ ^KbiKaLp + KpLa), ( 34 ) 


Finally 

Zap = Hap + iirq^eapp. (49) 

or 

Zap = ZiKpSa + ZllKaSp + ZiuLaSp + ZivLpSa (50) 

where 

Zj = {Hi — Airqn)] Zn = {Hi + 47rg7/); 

Zui = {H 2 - 477(7/); Ziv = {H 2 + 47rq/). (51) 

YaT[:iah\es:YTjji,KL,XT,i,ii,KL,Zijijijjv are the 
structure scalars of our distribution. 


Hap = Hi{SaKp + SpKa) + H2{SaLp + SpLa). ( 35 ) 
Let us now introduce the following tensors 
Yap = Ra.psV'V^ 


I 


^ 0/3 = ^Va.^‘^KppsV''V\ 


ZaB = ^eaepRsp"^V\ 


( 36 ) 

( 37 ) 

( 38 ) 


where R^p^sZZl^^P’^sRap^''■ 
Or, using 


Yap = -Yrhap + 2(2^^ + Yii){KaKp — -hap) 



+ 2^^^^ Yi){LaLp — —hap) 



+ YKL{KaLp + KpLa), 

(39) 

with 

YT = M^^ + 3P), 

(40) 


Yi=£i- d/rB/, 

(41) 


Yu = £ii — dTrll//, 

(42) 


Ykl = £kl - ‘^'^Hkl- 

(43) 

Xap 

= -^XThap + + Xu){KaKp — 

+ -^{‘^Xjj + Xi){LaLp — -hap) 

'^hap) 


+ XKL{KaLp + KpLa), 

(44) 

with 

Xt = SiTfi, 

(45) 


Xi = —£i — 4/7II/, 

(46) 


Xu = —£ll — 47711//, 

(47) 


Xkl = —£kl — ^tiHkl- 

(48) 


C. The super—Poynting vector 

An important role in our discussion is played by the 
super-poynting vector. Indeed, we recall that we define 
a state of intrinsic gravitational radiation (at any given 
point), to be one in which the super-Poynting vector does 
not vanish for any unit timelike vector [2ll - l23l |. Then 
since the vanishing of the magnetic part of the Weyl ten¬ 
sor implies the vanishing of the super-Poynting vector, it 
is clear that FRW does not produce gravitational radia¬ 
tion. It is also worth recalling that the tight link between 
the super-Poynting vector and the existence of a state of 
radiation, is firmly supported by the relationship between 
the former and the Bondi news function [2J (see [2^ 
for a discussion on this point). 

Then from the definition of the super-Poynting vector, 

Pa = eapj {Y^Z^^ - XJZ^P) , (52) 

we obtain 

Pa = PiKa + Pi I La, (53) 

with 

Pi = ^{^Yii + Yp- 2Xii - A/) + Hi{Ykl - Xkl) 

H-^ [2Tr + 2 At — Xj — Yj] — A'Kqii{XKL + Ykl), 

Pii = “^(22^3 + Xii — Yu — 2Yi) + H2 {Xkl — Ykl) 

— 4:TTqi{YKL + Xkl) 

+ [2yT + 2 At - Xu - Yu]. (54) 

Both components have terms not containing heat dis¬ 
sipative contributions. It is reasonable to associate these 
with gravitational radiation. Also, note that both com¬ 
ponents have contributions of both components of the 
heat flux vector. 

There is always a non-vanishing component of P^, on 
the plane orthogonal to a unit vector along which there is 
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a non-vanishing component of vorticity (the 6 — r- plane). 
Inversely, vanishes along the (/)-direction since there 
are no motions along this latter direction, because of the 
reflection symmetry. 

We can identify three different contributions in (l54l) . 
On the one hand we have contributions from the heat 
transport process. These are independent of the mag¬ 
netic part of the Weyl tensor, which explains why they 
remain in the spherically symmetric limit. 

On the other hand we have contributions from the 
magnetic part of the Weyl tensor. These are of two kinds: 
a) contributions associated with the propagation of grav¬ 
itational radiation within the fluid, b) contributions of 
the flow of super-energy associated with the vorticity on 
the plane orthogonal to the direction of propagation of 
the radiation. Both are intertwined, and it appears to be 
impossible to disentangle them through two independent 
scalars. 

As mentioned before, both components of the heat flux 
four-vector, appear in both components of the super- 
Poynting vector. Observe that this is achieved through 
the Xkl + difi terms in (iMl) , or using (l43l l48l) , through 
^kl- Thus, n^L couples the two components of the 
super-Poynting vector, with the two components of the 
heat flux vector. 


III. THE EQUATIONS 

We shall now deploy the whole set of equations for the 
variables defined so far. 


A. The heat transport equation 


We shall need a transport equation derived from a 
causal dissipative theory (e.g. the Miiller-Israel-Stewart 
second order phenomenological theory for dissipative flu¬ 
ids [li-lil). 

Indeed, the Maxwell-Fourier law for radiation flux 
leads to a parabolic equation (diffusion equation) which 
predicts propagation of perturbations with infinite speed 
(see [13-[13 and references therein). This simple fact 
is at the origin of the pathologies ISJ fo und in the ap¬ 
proaches of Eckart 341 and Landau [13 for relativistic 
dissipative processes. To overcome such difficulties, var¬ 
ious relativistic theories with non-vanishing relaxation 
times have been proposed in the past [26M29I . ■ The 
important point is that all these theories provide a heat 
transport equation which is not of Maxwell-Fourier type 
but of Cattaneo type [ 11 ] , leading thereby to a hyperbolic 
equation for the propagation of thermal perturbations. 

A fundamental parameter in these theories is the re¬ 
laxation time T of the corresponding dissipative process. 
This positive-definite quantity has a distinct physical 
meaning, namely the time taken by the system to return 
spontaneously to the steady state (whether of thermo¬ 
dynamic equilibrium or not) after it has been suddenly 


removed from it. Therefore, when studying transient 
regimes, i.e., the evolution between two steady-state sit¬ 
uations, T cannot be neglected. In fact, leaving aside that 
parabolic theories are necessarily non-causal, it is obvi¬ 
ous that whenever the time scale of the problem under 
consideration becomes of the order of (or smaller than) 
the relaxation time, the latter cannot be ignored, since 
neglecting the relaxation time amounts -in this situation- 
to disregarding the whole problem under consideration. 

Thus, the transport equation for the heat flux reads 

[lamiiii, 


rKq-pV^ + q^ = -{T,+Ta,)-\nT^ 

’“(55) 

where r, k, T denote the relaxation time, the thermal 
conductivity and the temperature, respectively. 
Contracting (|55|) with we obtain 


^ ^ ( GTt + A^Te , ^ 

-r \qii,t + Aqifl) + qii = -— , "b AT an 


A 


A VVA2B2r2 -hG2 
2 


(56) 


where (EH), has been used. 

On other hand, contracting (1551) with AT^, we find 



AqjjQ) + qj 


-^{Tr + BTai) 

KT^qijV^^ 

2 


(57) 


It is worth noting that the two equations above are cou¬ 
pled through the vorticity. This fact entails an interesting 
thermodynamic consequence. Indeed, let us assume that 
at some initial time (say t = 0) and before it, there is 
thermodynamic equilibrium in the 9 direction, this im¬ 
plies qii = 0, and also that the corresponding Tolman’s 
temperature [^ is constant, which in turns implies that 
the term within the round bracket in the first term on 
the right of (1551) vanishes. Then it follows at once from 
(l56l) that: 


qii,t = -Anqi, (58) 

implying that the propagation in time of the vanishing 
of the meridional flow, is subject to the vanishing of the 
vorticity and/or the vanishing of heat flow in the r- di¬ 
rection. 

Inversely, repeating the same argument for (1571) we ob¬ 
tain at the initial time when we assume thermodynamic 
equilibrium, 


qi^t = Aflqii. (59) 

Thus, it appears that the vanishing of the radial com¬ 
ponent of the heat flux vector at some initial time, will 
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propagate in time if only, the vorticity and/or the merid¬ 
ional heat flow vanish. 

In other words, time propagation of the thermal equi¬ 
librium condition, in either direction r or 9, is assured 
only in the absence of vorticity. Otherwise, it requires 
initial thermal equilibrium in both directions. 

This result is a clear reminiscence of the von Zeipel’s 
theorem (40j| . 


B. The equations for the metric functions, the 
kinematical variables and the Riemann tensor 
components. 


Let us hrst recall the decomposition of the covariant 
derivative of the four-velocity in terms of the kinematical 
variables given by: 

— ^a/9 “b f^a/3 “b ^^ 0 / 30 , (6^) 

which entails all the equations (fT^ . (fT^ . (120]) . ([24|) . 

Now, if we regard the above expression as a first order 
differential equation relating the kinematical variables 
with first order derivative of the metric functions, and 
look for its integrability conditions, we find 


Va-l3-u — Va-u-,P = 

From this last equation the following equations are ob¬ 
tained, by projecting with different combinations of the 
tetrad vectors: 

An evolution equation for the expansion scalar (the 
Raychaudhuri equation) 

+ ^0" + 2((t2 - - a% + + 3P) = 0 (62) 

O 


where 2cr^ = 

An equation for the evolution of the shear tensor: 


~ ^2(7 -f ^;6) ^ap “b ^a^P ^a^P 

- + Eap - dTTlIa^ = 0. 


(63) 


An equation for the evolution of the vorticity tensor: 


h^h'pQ,^v-,5y^ + —012(1/3 + 2(T^[(in^^j — — 0. 

(64) 

Two constraint equations relating the kinematical vari¬ 
ables and their derivatives with the heat flux vector and 
the magnetic part of the Weyl tensor: 




+ (o’(i/3 + ^ap) = Snqa, 


(65) 


C. The conservation equations 

The conservation law = 0, leads to the following 
equations: 


P'\ay°‘ + (m + P)y> + ^(20”/ + O’//)!!/ 

+ i(2(T//-b O’/)!!//-b -b g“a(i = 0, (67) 


(m + P)aa + (^P.p -I- -I- qp-^jV^ 


+ 1 -^Qhap + O'aP + ^aP ) <7^ — 0. 


( 68 ) 


D. The Bianchi identities 

Next, if we regard (EH) as a system of differential equa¬ 
tions of first order, relating the Riemann tensor compo¬ 
nents with the kinematical variables and their deriva¬ 
tives, and look for their integrability conditions, we are 
lead to the Bianchi idenitities, which together with (1281) . 
lead to the following set of equations: 

An evolution equation for the electric part of the Weyl 
tensor 

^{a^'p)R' l^‘'',^y^ + QEap + hapE^i,a^’^ — 3ii/^((iO’^^ 

47r : 

- 47r(// -b P)cTap - —0!!a/3 “ 47r7i'/^/i^)!!^(,;5 - 47rcr^((i!!j^^ 

47r 

- 47rn^(^n,g)^87ra((i(7^) + y -b a^q^ + g)/) Kp 

- (69) 

A constraint equation for the spatial derivatives of the 
electric part of the Weyl tensor 

+ SH^pco^ = 

- 47Thih>^^n^,,^ - Air 

A constraint equation for the spatial derivatives of the 
magnetic part of the Weyl tensor 

{^(JasE^P + S^aSEp'j + (EElyn ~ = 

Anip, + P)nape^°‘^ 

+ 47r [qa-,p + !!kc((c’’^^ -b p)\ ^k°‘^^ (71) 


ehi-aP^ + snj]qp, (70) 


- ,u 1 / ^ ^ uK-y/iir TT /r.r.\ Au ovolutiou equatlou for tho magnotic purt of tho Woyl 
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+ 

+ 

+ 


+ -QH^ 

K,5y^KKi - (^-5 + 

(a^s + n^s) 


47r 

Y 

47r 


/ + 47rn„^.;/3/i;^e 


QhS^af5 4“ ^ck(^k/3 4“ ^k/S 4“ 


(72) 


Equations (l55|) . (I62l) - (l72]) form the full set of equations 
for the variables of our problem. However, the following 
remarks are in order at this point: 

• Obviously, not all of these equations are indepen¬ 
dent, however depending on the problem under con¬ 
sideration, it may be more advantageous to use one 
subset instead of the other, and therefore here we 
present them all. 

• The scalar equations obtained by projecting the 
above system, on all possible combinations of 
tetrad vectors, are deployed in the Appendix B of 

la 

• The obtained equations are of first order, unlike 
the Einstein equations, which are differential equa¬ 
tions of second order for the metric functions. This 
reduction is achieved by enlarging the number of 
variables and equations. 

• In the case of specific modeling, another important 
question arises, namely: what additional informa¬ 
tion is required to close the system of equations? 
It is clear that information about local physical as¬ 
pects of the source (e.g. equations of state and/or 
information about energy production) are not in¬ 
cluded in the set of deployed equations and there¬ 
fore should be given, in order that metric and mat¬ 
ter functions could be solved for in terms of initial 
data. 


IV. THE EFFECTIVE INERTIAL MASS 
DENSITY OF THE DISSIPATIVE FLUID 

In classical dynamics the inertial mass is defined as the 
factor of proportionality between the three-force applied 
to a particle (a fluid element) and the resulting three- 
acceleration, according to Newton’s second law. 

In relativistic dynamics a similar relation only holds 
(in general) in the instantaneous rest frame (i.r.f.), since 
the three-acceleration and the force that causes it are not 
(in general) paralell, except in the i.r.f. (see for example 
111 ). 

We shall derive below, an expression for the effective 
inertial mass density for our dissipative fluid distribution. 


By “effective inertial mass” (e.i.m.) density we mean 
the factor of proportionality between the applied three- 
force density and the resulting proper acceleration (i.e., 
the three-acceleration measured in the i.r.f.). 

As we shall see, the obtained expression for the e.i.m. 
density contains a contribution from dissipative vari¬ 
ables which reduces its value with respect to the non- 
dissipative situation. Such decreasing of e.i.m. den¬ 
sity was brought out for the first time in the spherically 
symmetric self-gravitating case in [d^. Afterwards this 
effect was also detected in the axially symmetric self- 
gravitating case (d^ . for slowly rotating self-gravitating 
systems Iddl . and under other many different circum¬ 
stances [dhl - l^ . 

It is perhaps worth noticing that the concept of effec¬ 
tive inertial mass is familiar in other branches of physics, 
thus for example the e.i.m. of an electron moving under 
a given force through a crystal, differs from the value cor¬ 
responding to an electron moving under the same force in 
free space, and may even become negative (see 1^1. 

Combining the equations (IB51) and (1551) we obtain 


ifi + P) (I - a) a„ = 
{o'a/3 4- fla/3)q^ 4-VqT 






0 ® f Qa, 


(73) 


an expression which takes the desired, ’’Newtonian”, 
form. 


Force=e.i.m. x acceleration(proper), 

where Dtf = and a = • 

The factor multiplying the four acceleration vector rep¬ 
resents the effective inertial mass density. Thus, the ob¬ 
tained expression for the e.i.m. density contains a contri¬ 
bution from dissipative variables, which reduces its value 
with respect to the non-dissipative situation. 

From the equivalence principle it follows that the “pas¬ 
sive” gravitational mass density should be reduced too, 
by the same factor. This in turn might lead, in some 
critical cases when such diminishing is significative, to a 
bouncing of the collapsing object. 

It should be observed that causality and stability con¬ 
ditions hindering the system to attain the condition 
a = 1, are obtained on the basis of a linear approxima¬ 
tion, whose validity, close to the critical point {a = I), is 
questionable [s^- 

At any rate, examples of fluids attaining the critical 
point and exhibiting reasonble physical properties have 
been presented elsewhere [13, HH . 

In order to evaluate a, let us turn back to c.g.s. units. 
Then, assuming for simplicity fi + p 2/i, we obtain 


_ M[r] , 0-42 
r(fJ-+p) ^ 


( 74 ) 
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where [k], [T], [t], [fj] denote the numerical values of 
these quantities in erg. s~^ cm~^ K, s and g. cm~^, 
respectively. 

Obviously, this will be a very small quantity (compared 
to 1), unless conditions for extremely high values of k and 
T are attained. 

At present we may speculate that a may increase sub¬ 
stantially (for a non-negligible values of r) in a pre¬ 
supernovae event 

Indeed, at the last stages of massive star evolution, 
the decreasing of the opacity of the fluid, from very high 
values preventin g th e propagation of photons and neu¬ 
trinos (trapping [56|), to smaller values, gives rise to ra¬ 
diative heat conduction. Under these conditions both k 
and T could be sufficiently large as to imply a substan¬ 
tial increase of a. Indeed, the values suggested in 
([k] ft! 1037; [T] ft 1013; [t] ft lO-'i; [g] ft lO^^ ) lead to 
a ft 1. The obvious consequence of which would be to 
enhance the efficiency of whatever expansion mechanism, 
of the central core, at place, because of the decreasing of 
its e.i.m. density. At this point it is worth noticing that 
the relevance of relaxational effects on gravitational col¬ 
lapse has been often exhibited and stressed (see [ssl - l^ . 
and references therein) 

It is also worth noticing that the inflationary equation 
of state (in the perfect fluid case) fi + P — 0, is, as far as 
the equation of motion is concerned, equivalent to a = 1 
in the dissipative case (both imply the vanishing of the 
e.i.m. density). 

Finally, it is worth stressing that it is the first term 
on the left and the second on the right, in (1551) the di¬ 
rect responsible for the decreasing in the e.i.m density. 
Therefore any hyperbolic dissipative theory yielding a 
Cattaneo-type equation in the non-relativistic limit, is 
expected to give a result similar to the one obtained here. 

V. SOME PARTICULAR CASES 

In what follows we shall consider some particular cases, 
where some variables (e. g. the shear) are assumed to 
vanish. We do so, on the one hand for simplicity, and 
on the other, in order to bring out the role of some spe¬ 
cific variables. However, it should be kept in mind that 
such kinds of “suppressions” may lead to inconsistencies 
in the set of equations. This is for example the case of 
“silent” universes [H, , where dust sources have van¬ 

ishing magnetic Weyl tensor, and lead to a system of 1+3 
constraint equations that do not seem to be integrable in 
general (^ . In other words for any specific modeling, 
the possible occurrence of these types of inconsistencies 
should be carefully considered. 

A. The shear free case 

This case has been analyzed in detail in 0 . Below we 
summarize the main results obtained under the shear- 


free condition. 

• For a general dissipative and anisotropic (shear 
free) fluid, vanishing vorticity, is a necessary and 
sufficient condition for the magnetic part of the 
Weyl tensor to vanish. 

• Vorticity should necessarily appear if the system 
radiates gravitationally. This further reinforces the 
well established link between radiation and vortic¬ 
ity. 

• In the geodesic (shear-free) case, the vorticity van¬ 
ishes (and thereof the magnetic part of the Weyl 
tensor). No gravitational radiation is produced. A 
similar result is obtained for the cylindrically sym¬ 
metric case, suggesting a link between the shear of 
the source and the generation of gravitational radi¬ 
ation. 

• In the geodesic (non-dissipative) case, the models 
do not need to be FRW, however the system relaxes 
to the FRW spacetime (if 0 > 0). Such tendency 
does not appear for dissipative fluids. 


B. The perfect, geodesic fluid 

In Q we have considered the case of perfect and ge- 
doesic fluid, without assuming ab initio the shear-free 
condition. As the result of such study we have found 
that: 

• All possible models compatible with the line ele¬ 
ment © and a perfect fluid, are FRW, and accord¬ 
ingly non-radiating (gravitationally). Both, the 
geodesic and the non-dissipative, conditions, are 
quite restrictive, when looking for a source of grav¬ 
itational waves. 

• Not only in the case of dust, but also in the ab¬ 
sence of dissipation in a perfect fluid, the system is 
not expected to radiate (gravitationally) due to the 
reversibility of the equation of state. Indeed, radi¬ 
ation is an irreversible process, this fact emerges 
at once if absorption is taken into account and/or 
Sommerfeld type conditions, which eliminate in¬ 
ward traveling waves, are imposed. Therefore, the 
irreversibility of the process of emission of gravi¬ 
tational waves, must be reflected in the equation 
of state through an entropy increasing (dissipative) 
factor. 

• Geodesic fluids not belonging to the class consid¬ 
ered here (Szekeres) have also been shown not to 
produce gravitational radiation. This strengthens 
further the case of the non-radiative character of 
pure dust distributions. 


C. The dissipative, geodesic fluid 

From the results discussed above, it becomes clear that 
the simplest fluid distribution which we might expect to 
be compatible with a gravitational radiation, is a dissipa¬ 
tive dust under the geodesic condition. Such a case was 
analyzed in 

The two possible subcases were considered separately, 
namely: the fluid distribution is assumed, from the be¬ 
ginning, to be vorticity-free, or not. 

In the former case, it is shown that the vanishing 
vorticity implies the vanishing of the heat flux vector, 
and therefore, as shown in Q, the resulting spacetime is 
FRW. 

In the latter case, it is shown that the enforcement 
of the regularity conditions at the center, implies the 
vanishing of the dissipative flux, leading also to a FRW 
spacetime. 

Thus all possible models, sourced by a dissipative 
geodesic dust fluid, belonging to the family of the line 
element considered here, do not radiate gravitational 
waves during their evolution, unless regularity conditions 
at the center of the distribution are relaxed. There¬ 
fore physically acceptable models require the inclusion 
of, both, dissipative and anisotropic stresses terms, i.e. 
the geodesic condition must be abandoned. In this case, 
purely analytical methods are unlikely to be sufficient to 


arrive at a full description of the source, and one has to 
resort to numerical methods. 


VI. OPEN ISSUES 

Below we display a partial list of problems which we 
believe deserve some attention: 

• How could one describe the “cracking” (splitting) of 
the configurations, in the context of this formalism 
? 

• We do not have an exact solution (written down 
in closed analytical form) describing gravitational 
radiation in vacuum, from bounded sources. Ac¬ 
cordingly, any specific modeling of a source, and 
its matching to an exterior, should be done numer¬ 
ically. 

• It should be useful to introduce the concept of the 
mass function, similar to the one existing in the 
spherically symmetric case. This could be relevant, 
in particular, in the matching of the source to a 
specific exterior. 

• What is the behaviour of the system in the quasi¬ 
static approximation? Would be there gravita¬ 
tional radiation in this case? 
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